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Abstract 

We consider the two-phase Stefan problem ut = a(u) where a(u) = u + 1 for 
u < —1, a(u) = for —1 < u < 1, and a(u) = u — 1 for u > 1. We show that if u 
is an L 2 distributional solution then a(u) is continuous. 

1 Introduction 

In this paper we will discuss the regularity of weak solutions to the two-phase Stefan 
problem 

| = A«(«) (1.1) 

in a domain Q C M n x (0, T), for some T > 0. Here a(u) — if — 1 < u < 1, = u — 1 
for it > 1, and = m + 1 for u < — 1. 

We will show that if u G L^ oc (f2) is a solution in the sense of distributions of (jl.lj) 
(denned precisely below) then a{u) is continuous. In the case when u > 0, u is a solution 
of the one-phase Stefan problem and Andreucci and Korten |AnKoj (see also Korten 
|Kop have shown that if u G Lj oc , then a(u) is continuous. Although we believe this to 
be true in the two-phase case, we have not been able to obtain this generality and must 
assume u G Lf oc . 

Under the assumption that u is bounded and Va(u) G L 2 , Caffarelli and Evans 
CaEj showed that a(u) is continuous. Similar results for more general singular parabolic 
equations were shown by Sacks [S], Ziemer [Z] and by DiBenedetto jDiBj . We will assume 
these results. We will show that a locally L 2 weak solution of (jl.lj) satisfies the hypotheses 
of these results (of any of these authors) to conclude the continuity of a(u). 

Related to this equation is the porous medium equation u t = Au m , m > 1. This has 
been studied extensively by many authors, but we mention in particular the regularity 
result of Dahlberg and Kenig jDKj who showed that a nonnegative solution to the 
porous medium equation is a.e. equal to a continuous function. The methods in this 
present paper are descendants (via the work of Andreucci and Korten) of the methods 
of Dahlberg and Kenig found in jDKj . However, the fact that we are working with 
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solutions which can be both positive and negative complicates matters. To achieve our 
results we will perform numerous integrations by parts and cannot determine the sign 
of the resulting boundary terms as in the one-phase case. Consequently we devise a 
different strategy and introduce new ideas and techniques. 

Equation (jl.ljl is a formulation of the two-phase Stefan problem, describing the flow 
of heat within a substance which can be in a liquid phase or a solid phase, and for which 
there is a latent heat to initiate phase change. This allows for the presence of a "mushy 
zone", that is, a region which is between the liquid and solid phases. In this model u 
represents the enthalpy and a(u) the temperature. We have assumed that the thermal 
conductivity in both the solid and liquid phases is the same. These conductivities are 
determined by the slope of the function a(u) in the regions u > 1, and u < — 1. The 
results below all continue to hold (with minor modifications) if the slope of a(u) differs 
in these regions. 

We now state our main result. Suppose u G L 2 oc (Q) where Q is a domain contained 
in R 71 x (0, T). We consider distributional solutions of the equation u t = Aa(u), that is, 
u which satisfy 

a(u) Aip + mptdx dt = 



n 

for every ip G C°° with compact support in Q. 

Theorem 1.1. Suppose u G Lf oc (Q) is a solution of u t = Aa(u). Then a(u) is a.e. equal 
to a continuous function. 

We do not expect, in general, such a result for u. As noted in Korten |Kolj . the 
solution to the Cauchy problem ut = Aa(u) on with initial data < ui(x) < 1 is 
just u(x, t) = Ui(x). Thus, we cannot expect u(x, t) to be any smoother than ui(x). 

The paper is structured as follows. In section 2 we prove energy estimates for weak 
solutions of the two phase problem. These show that Va(u) and a(u) t exist locally in L 2 . 
In section 3, we show that |q;(u)| is subcaloric. An immediate consequence is that a(u) 
is locally bounded. This, combined with the energy estimates and previously mentioned 
theorem of DiBenedetto |DiB| (or others mentioned above) gives the continuity of a(u). 



Throughout, the letter C will denote a constant which may vary from line to line. 
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The second author would like to thank the National University of Ireland at Galway for 
their hospitality during part of this work. 

2 Energy Estimates 

We establish that a(u) has derivatives which are locally in L 2 . 

Theorem 2.1. Suppose Q C 1R™ +1 and u G L 2 ^ (tt) is a distributional solution of 
u t = Aa(u) on Q. Suppose r < R, T < t < ti < T\, set uj = (to?^i) x B(x ,r) and 
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Co = (T ,Ti) x B(xq,R), and suppose the closure of Co is contained in Q. Then Va(«), 
a(u) t exist in L 2 (u) and there exists a constant C , depending only on uo and Co such that 



Va{u)\ 2 dxdt <C u 2 dxdt (2.1) 



and 



8 M 2 



dxdt<C // u 2 dxdt (2.2) 



Proof. Let ip m (y,s) = p m (y)T m (s), m = 1,2,... where p m , r m are smooth mollifiers, 
radial, centered at 0, compactly supported, and tending to <5o- For (x, t) G Q and m 
sufficiently large (depending on (x, t)), ip m (x — y,t — s) is a test function supported in Q 
and thus 

u(y, s)-^(x -y,t-s) + a(u(y, s))Ap m (x -y,t- s)dy ds = 0. 

In the course of the proof, we will need to define three nested domains between uo and 
Co. To simplify notation, set u)\ = oo, oo 5 = Co and we will define co>2, uj 3 and oo A with 
Wi C w 2 C W3 C W4 C uj 5 . 
For m = 1, 2, 3, . . . set 

u m {x,t) = J J u(y, s)xu> s {y, s)(p m {x - y,t - s)dyds. 

Then for all (x,t) and m, \u m (x,t)\ < M(uXui r A x ^)) where M denotes the Hardy- 
Littlewood maximal function (in both the variables (x,t)). 

Choose a, < a < and b, <b< *i±pL suc h that 

\M(u X ^)(x,a)\ 2 dx<C [[ \M{u X ^){x,t)\ 2 dxdt. (2.3) 

B(x ,R) J Jus 



with a similar inequality for b. 

In a similar fashion, set w m = cx(u)xu s *ipm- Define uo$ = B (x , ^-^f^) x ( 
Then on to 4., jr,u m — Aw m = for all m sufficiently large. Using cylindrical coordinates 



3Tb +f ti+3Tj ' 
4 ' 4 



we can choose an r\, < ri < r+ 4 3R so that 

\M{a{u) X ^)\ 2 da<C ff 

zo,n)x(T ,Ti) J Jui 

Then by ()2.3|) . for all sufficiently large m 



\M(a(u)xu; 5 )(x,t)\ 2 dxdt. (2.4) 

dB(x ,n)x(T ,Ti) J Jui 5 



[ „ v , ...... . , 

JB(xo,n) JB(x ,R) 



u m (x,a) 2 dx< / |M(«x £Js )(a;, a)| 2 <ix 



<C// |M(MXc 5 )(x,t)| 2 rfxrft (2.5) 
<C I I u 2 dx dt 
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Likewise, by (|2~1) 



w 2 m da< / \M(a(u) X ^)\ 2 da 

dB(x ,ri ) x (a,6) JaB(x ,n)x (T ,Ti ) 



<C// |M(a(u)x U6 )| 2 da;cft (2.6) 
< C // H«)| 2 dxdt. 

Let at m (s) be a smooth regularization of a(s) such that a m (s) = a(s) for \s\ > 1 + — , 
a m (s) is strictly increasing and a m (s) ^ except for s — 0. Put co> 3 = I?(a;o,Ti) x (a, 6). 
Let f m be a solution to 

v t = Aa m (v) on u; 3 

v (x, a) = u m (x, a) x £ B{xQ,r\) 
a m (v) =w m on dB(x Q ,r 1 ) x (a,b). 

Choose (p(x) so that (f> = on dB(x , ri), A</> = 1 on B(xq, ri). Then < on B(xq, ri) 
and |^ = ci > on dB(xo,ri), where n is the outward normal and ci is a constant 
depending only on ri and the dimension. 
By Green's theorem we have 

a m (v m )) 2 A<f)dx 

B(x ,ri) 

A(a m (v m )) 2 (f)dx + / (a m (u m )) 2 —da - / <j)—[a m (v m )} 2 da 

B(x ,ri) JdB(x ,ri) O n J dB{x ,r-i) c,n 

= 2 Aa m (v m )a m (v m )(j)dx + 2 | Va m (v m )\ 2 (j)dx + ci / a m (v m ) 2 da 

J B(xo,ri) J B(xQ,ri) JdB(xo,r\) 

<2 v mt a m (v m )(f)dx + ci / a m (w m ) 2 <io- 

J B(xo,r\) JdB(xo,ri) 

= 2 4 / A m (w m )0rfx + ci / a m (v m ) 2 da 

dt JB(x ,ri) JdB(x ,ri) 

where A m is an antiderivative of a m . Integrate from a to 6 to obtain 



na m (v m ) 2 A(j)dxdt < 2 / A m (t; m (x, 6))0c£r - 2 / 
i(a;o,ri) J B(x ,n) J B{ 



A m (v m (x,a))(f)dx 

I B(x ,r 1 ) J B{xQ,r{) J B(x ,r 1 ) 

+ / a m (v m ) 2 dadt. 

J a J dB(xo,ri) 
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Now < A m (x) < x 2 , cf) < 0, A0 = 1 and recalling u 3 = B(x , T\) x (a, b), this yields: 
a m (v m ) 2 dx dt < C / v m (x, a) 2 dx+ / / a m (v m ) 2 dadt 

U13 J B(xQ,n) J a J dB(xQ,r\) 

C / u m (x, a) 2 dx+ / / w 2 dadt (2.7) 



-B(xo,ri) ./a J dB(xo,ri) 



<C u 2 dxdt 



where for the last inequality we have used ()2.5|) and ([2.6)1 . Let ?/>(x) be a nonnegative 
C °°(M n ) function such that ip = 1 on B(ar , i[> = outside B(a? , ^Ip). To simplify 

notation set B(x , ^±£) = B x , B(z , ^) = S 2 . Then 

ipa m (v m )v mt dx = / V ; «m(f m )Aa m (f m )(ix 

B 2 ^^2 



Vi>-Va m (v m )a m (v m )dx- / ^|Va m (w m )| cte 
b 2 Jb 2 

Aipa m (v m ) 2 dx - / ^|Va m (t; m )| 2 cfo. 

B 2 JB 2 



Rearrange and integrate from a to h to obtain 

il){Va m {v m )) 2 dxdt = ]- [ [ Aip\a m (v m ) 2 \dxdt- f f i>—A m {v m )dxdt 

a J B 2 * J a J B 2 J a J B 2 dt 

< C / / a m {v m ) 2 dxdt + / ^(^An-Kn,^, a)) da; 
<C / / a m (v m ) 2 dxdt + / v m (x,a) 2 dx 

Ja J B 2 J B 2 

<C a m (v m ) 2 dxdt + / v m (x,a) 2 dx 

J J U13 J B(xa,r\) 

<C u 2 dx dt 

J J w$ 

(2.8) 

where we have used ()2.7|1 and (|2.5jl and the definition of v m for the last inequality. 

We now seek a similar estimate for the t derivative. Let rj(x) be a nonnegative C™(M. n ) 
function such that r] = 1 on B(x ,r), r) = outside £>i and so that ||^=||oo < oo. Note 
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that a m (v m ) t = a' m (v m )v mt and < a' m < 1 so that {a m (v m ) t ) 2 < a m (v m ) t v mt . Then 
rj(a m (v m )t) 2 dx < / ■qa m (v m ) t v nH dx 

Bi JBx 

= / T]a m {v m ) t ^a m {v m )dx 

= - Vr? ■ Va m (v m )a m (v m ) t dx - / r]Va m (v m ) t ■ Va m (v m )d 
Jbj j b x 

= -[ Vr? • Va m {v m )a m {v m ) t dx - ]-^- [ r]\Va m (v m )\ 2 dx 

J B 1 * JBi 

Integrate from c to d, where c and d are to be chosen momentarily. We obtain 
r d r 

r)(a (v m ) t ) 2 dx dt 



x 




c J B 1 



< 



< 



/ / Vrj—pV 'a m (v m )a m (v m ) t dx 

Jc JBi 

\^a m (v m )\ 2 dx dt 



Bi Vv 



1 



dt+- r](x)\Va m (v m )(x,c)\ 2 dx 



lB x 




oo ^" c J Bi 

+ o [ v( x )\^ a m(v m )(x,c)\ 2 dx 

1 JBx 



r]{a m (v m ) t ) 2 dxdt I 

c JBi J 




Choose c m (depending on m), T °+ 3t ° < c m < to, so that 



(2.9) 



/ r](x)\Va m (v m )(x,c m )\ 2 dx < C / ijj\Va m (v m )\ 2 dx dt. 

>Bi J a JB 2 



(2.10) 



Put d — ti, c — Cm in ()2.9j) . Then recalling that ip = 1 on B\, and using (j2.1()j) and ()2.8j) 

we have 




Cm </i?l 

V?7 



< 



< 



r](a m (v m ) t ) 2 dxdt 




vfj\V a m (v m )\ 2 dx dt 



X) \ jc m J B2 

2 




r}(a m {v m )t) 2 dx dt J +C u 2 dx 

I Cm, J B\ J J J U)$ 



CO 




r](a m (v m )t) 2 dx dt I +C u 2 dxdt 

I C m J Bl / J J 



from which it follows that 

'■ti 



nr](a m (v m ) t ) 2 dx dt < C // u 2 dxdt 
. . — h J J ujf* 



and consequently 

{a m {v m )t) 2 dxdt < C II u 2 dxdt. (2-H) 



Thus, recalling uj = u>±, cu = u>5, ()2.8|) and ()2.11|) give 

Va m (v m )\ 2 dx dt < C / / u 2 dxdt and / / (a m (v m ) t ) 2 dx dt < C / / u 2 dxdt. 



(2.12) 

To obtain (j2.1J) and ()2.2|) we will need to take limits. We first remark that with more 
care, similar estimates could be obtained with any compact set K d u 3 replacing u = uj\ 
on the left hand side of the inequalities in (|2.12j) ; naturally, the constants on the right 
hand side depend on the position of K within u 3 . Thus, from (J2.7)) and this observation, 
we have: 

/ / a m (v m ) 2 dx dt < C u 2 dxdt, \ \ \Va m {v m )\ 2 dx dt < C(K) / / u 2 dxdt 

J J W3 J J to J J K J J LO 

and jj -^a m (v m ) dx dt < C(K) JJ u 2 dx dt 

(2.13) 

for every compact K C u 3 . 

By Rellich-Kondrachov there exists a subsequence {ot mk {v mk )} of {a m {v m )} (which 
we still write as {a m (v m )}) and h e L 2 (u 3 ) such that a m {v m ) — > h in L 2 {K) for every 
compact set K C uj%. By taking subsequences, if necessary, we also may assume this 
convergence is a.e. By weak compactness, and again, by taking subsequences, we may 
assume that a m (v m ) — ► h weakly in L 2 (u; 3 ). Equation ()2.13|) implies that the L 2 (co 3 ) 
norms of the v m are uniformly bounded, hence there exists a subsequence, (still denoted 
by v m ) such that v m — > v G L 2 (u 3 ) weakly. 

We claim that a(v) = h. First note that \\a m — — > 0, so that for a.e x G u 3 , 
(%(v m ) h. Consider the set where h > 0. Then for a.e x in this set, a(v m (x)) — > 
h(x) > 0, and hence v m (x) — > h(x) + 1. Thus, v(x) = h(x) + 1 for a.e x in the set 
where h(x) > 0. Similarly, on h < —1, v(x) = h(x) — 1 a.e. On the set h(x) = we 
must have —1 < liminfv m (x) < limsupf m (x) < 1 a.e. To see this consider an x at 
which there exists a subsequence v mk (x) which converges to yo ^ [—1,1]. Then for this 
x, a(v mk (x)) — > a(yo) 7^ which implies a(v m (x)) h(x). Thus, —1 < liminft> m (x) < 
limsupf m (x) < 1 a.e. on h(x) = 0, and hence —1 < v(x) < 1 a.e. on h = 0. We 
conclude that a(v) = h a.e. 

Summarizing, we have v m — > v weakly in L 2 (u 3 ) and a m (v m ) — > a(v) weakly in 
L 2 (u 3 ), a.e. on u 3 and in L 2 {K) for every compact subset K of oj 3 . To finish the proof 
we show that a(u) = a(v) a.e. on u 3 . Using integration by parts, and recalling that 
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a m( v m) — w m on dB(x , r x ) x (a, b) , we compute 



(v m - u m )(a m (v m ) - w m )dxdt 



a/ v mt (x,r) - u mt {x,T)dT(a m {v m {x,t)) - w m (x,t))dx dt 
!(xo,ri) J a 

n\ V(a m (f m ) - w m )(x,r)drV(a m (v m ) -w m )dxdt 
-J(xo,ri) J a 




a J B(xo,ri) 



d_ 

It 



V(a m (w m ) - w m )dr 



dx dt 



V(a m (f m ) -w m )dr 



dx < 



'B(x ,n) 

We need to take limits as m —>■ oo in this inequality. Write 



(2.14) 



{v m - u m )(a m (v m ) - w m )dxdt = v m a m (v m )dx dt + v m (-w m )dx dt 

J J U3 J J LOZ 

—u m )(a m (v m ))dx dt + u m w m dx dt = I + II + III + IV 



^'3 



Since 

Urn ^ U 8b. e. and in L 2 (co>3), w m — > a(u) a.e. and in L 2 {ojz), v m — > v weakly, and 
&m(vm) — ¥ cx(v) weakly, we conclude 



II 



v(—a(u))dx dt, III 



-u)(a(v))dx dt, 



U>3 



and IV 



ua(u)dx dt. 



Expand out 



{a m (v m ) - a m (v))(v m - v)dx > 0, 



0J3 



take m — > oo (make use of the fact that \\a m — a Hoc — > 0) to conclude 

liminf / / a m (v m )v m dx > / / a(f)w(ix. 



"-'3 



This combined with the estimates for II-IV and (J2.14|) yields 



(v — u)(a(v) — a(u))dxdt < 0. 



Since the integrand of this is nonnegative, we conclude a(u) = a(v) a.e. on U3. This 
completes the proof of the theorem. □ 
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3 \oi(u)\ is subcaloric 

Theorem 3.1. \ot(u)\ is weakly subcaloric, that is, it satisfies 

-V\a(u)\Vr) + \a(u)\r) t dx dt > 



for any nonnegative 7] G W^' 2 (Q) 



Proof. Let < rj G W 1,2 (fi). For h > set 



if x > h 



<t>h{x) = < 



\x - 1 |<x</i 

ft 2 — 

if |x| < | 

fx + 1 -/i<x<-§ 

-1 if x < -h. 



Then i](j) h (a(u)) is supported in ||a(-u)| > |} and thus, ff(a(u) — u)[r}<fih(a{u))\ t dx dt = 
0. 

Then 

= yy a(u)[770 fe (a(u))] t - Va(u)V[r)(j> h ((x(u))]dxdt 

a(u)r) t <fih(oi{u))dx dt + a(u)r](f)' h (a(u))a(u) t dx dt 
J J Va(u)Vri(i) h (a(u))dxdt — J J V a(u)r](j)' h {a{u))V a{u)dx dt 



(3.1) 



I + 11 + III + IV 



We investigate each of these as h — > 0. As /i — > 0, 0ft(a('u)) — > sgn(o;(w)) so that 
/ — > JJ |a(u)|ntda;(it. To estimate 77, first note that 



Then 



and consequently, 



= ^X{| <H< ft}- 
/J = JJ a ^ v l x ^< a ( u )\< h }§i a ^ dxdt 



\II\ < 



Hu)\rj- X{ ^ <H u)\<h} 



d_ 

dt 



a(u) 



dxdt<2 / / VX { !i <a{ u)<h} 



d_ 

dt 



a{u)) 



dx dt. 



Since §ra(u) G L 2 loc {ti), II -> as /i -> 



To estimate note that when |a(w)| > ft, Vet («(«)) = V|«(m)|. And when 
§ < \a(u)\ < ft, 



\Va{u)(t) h {a{u))\ < \Va(u)\ [ -r\a(u)\ + 1 



Consequently, 



{\<\a{u)\<h} 



V 'a(w) V 'rj(p h {a{u))dx dt 



< I I \Va(u) 
'{|<Ku)|<fc} 



\Vi]\dx dt 



-» as h — »• since Va(it) G L 2 oc (fi)). 
Therefore, as ft — > 0, III —>■—]] V|a(w) \ Vrjdx dt. Note that we can write IV as 

2 



IV = JJ l Va ^' 2r? flh( a ( u )) dx dt = jj \^ a ( u )^V^X{h <]a(u)l<h ydxdt 
Thus, letting ft -> in (JSHJ) yields: 

= J J \a{u)\rj t dx dt — I f V|a(w)| Vrjdx dt — lim yy | Vct(u)| 2 — X{fc<[ a ( ti )[<h} 7 7 £ ^ c ^ 
from which the theorem follows. 



□ 



Remark. Suppose that instead of <ph as defined above, we defined 

1 



if x > h 
j-x — 1 if | < x < h 

h 2 

otherwise. 



Then following the computations as in (|3.1|) we obtain (|3.1|) with this version of <ph- 
In this case J — > JJ a(u) + f] t dx dt and as in the above case, // — > 0, and III — > 
- JJ Va(u) + Vr]dxdt. 
Now we may write 

/V = // |Va(M)+|2?? ^ X {f «*(«)<fc} dxc ^ 

We obtain: 

= a(u) + i] t dxdt — JJ Va{u) + Vrjdx dt — lim y^ |Va(w) + | 2 — X{A< Q! ( u ) < / l } 7 7^ a; ^ 
Thus, + is subcaloric. In a similar fashion, we may use the function 



\x + 1 —ft, < x < — | 



</>h(a;) 



if x < —h 
otherwise 
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and computations such as those above to obtain 



= 



// 



a(u) rjtdxdt 




Va(ii) Vrjdxdt + lim 

h— »0 



JJ\V*{u)-f- X{ 



h<a{u )<-±}Vdxdt 



to conclude that is supercaloric. 
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